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1 . Introduction 


Although the effect of heavy rain on aircraft performance was discussed 
as early as 1941 by Rhode (1), no serious studies of the relationship of heavy 
rain to aircraft safety were made until Luers suggested in 1981 that the tor- 
rential rain that often occurs at the time of severe wind shear might 
substantially increase the danger to aircraft operating at slow speeds and 
high lift in the vicinity of airports. While Luer’s ideas were not published 
until early 1983 (2), appropriate measures were taken by NASA to study the 
effect of heavy rain on the lift of wings typical of commercial aircraft. 

These tests, reported by Dunham, Bezos, Gentry, and Melson (3), were the sub- 
ject of a number of discussions between the senior author of this report and 
Nr. Earl Dunham of NASA during the fall of 1984. One of the aspects of these 
tests that seemed confirmed by the data was the existence of a "velocity 
effect" on the lift data. The data seemed to indicate that when all the nor- 
mal non-dimensional aerodynamic parameters were used to sort out the data, the 
effect of velocity was not accounted for, as it usually is, by the effect of 
dynamic pressure. Indeed, the measured lift coefficients at high lift indi- 
cated a drop'-off in lift coefficient for the same free-stream water content as 
velocity was increased. 

This behavior was explained at the time by the authors as being due to 
the variations in momentum deposition of the droplets that are splashed back 
into the air stream after shattering upon striking the airfoil surface. Since 
the higher the speed the smaller are the splashed back droplets, it follows 
that the higher the speed the closer to the airfoil will be the layer where 
these droplets are reaccelerated by the air stream. It was suggested that if 
this splashed-back momentum effect was large enough and took place close 
enough to the airfoil, the airfoil could stall. This would be a very serious 
condition if, indeed, it could occur. These ideas were put forward at a meet- 
ing of contractors working the heavy rain problem for NASA in April of 1985 


Subsequent to this meeting, funds were obtained from NASA to modify the 
A.R.A.P. ARB code to allow the computation of the effect of splashed back 
droplets on the location of separation on airfoils at high angle of attack to 
see if stall could be induced and under what conditions of speed and scale 
such a phenomenon might be hazardous to aircraft operation. A secondary ob- 
jective of this calculation was to obtain such information as could be gleaned 
to aid NASA in conducting properly scaled tests of airfoils in heavy rain. 

Bilanin (5) has given a more detailed review than is attempted here of 
the mechanisms that may be involved when aircraft operate in heavy rain. 


2 . Ejecta Scaling 


To derive an expression for the extent of the momentum-defect layer that 
results from droplet splash-back, we first determine the order of magnitude of 
the diameter of the droplets that splash back from a surface after drop im- 
pingement . 

We consider a drop of radius r Q impinging normally with velocity V on a 
frictionless surface unwetted by the drop. As the drop hits the surface it is 
deformed into a pancake whose instantaneous radius is a and whose thickness is 
t as shown in Figure 2.1. From continuity, we have the result that 


3 nr 0 = " a2t 


( 1 ) 


At any given time, the surface energy of the distorted pancake is 


o w (2tra 2 + 2nat) 


( 2 ) 


where o is the surface tension of water in contact with air. 


In view of (1), this may be written 
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Generally when the drop breaks up into droplets, t«a, so we write the surface 
energy approximately as 


~ nr 2 • ^ a 
3 0 t w 


(4) 
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Figure 2.1 Geometry of drop impingement. 


Now this surface energy cannot increase without bound as t becomes smaller, 
for the only energy available to supply this surface energy is the kinetic 
energy of the original drop. In an actual impingement process, energy is lost 
to both internal and external friction, but if we neglect these and other 
losses and say that droplets of diameter d^ = t must form where a portion of 
f of the original kinetic energy has been converted to surface energy, we may 
write 
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p V 2 

w 
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( 5 ) 


or 
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p V 2 
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( 6 ) 


The above expression gives the order of magnitude of droplet size that can be 
formed. From this expression we learn that the faster one flies through the 
air the smaller are the drops that are splashed back from the surface of an 
airfoil. Also, we see that the droplet splash-back size is expected to be 
relatively independent of the initial drop size r . 

Having obtained an expression for the size drop one expects to see 
splashed back, we may now derive an expression for a length that is typical of 
the distance over which the splashed drops are accelerated. Consider the 
deceleration of a particle by its drag: 
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where m and A are the mass and cross-sectional area of the particle, C_. is 
P P D 

the coefficient of drag, p is the air density, V is the magnitude of the 

El 

velocity, and t is time. We may write this, since ds = Vdt, as 


dV 

V 




( 8 ) 


There is clearly a characteristic length associated with the acceleration 
process, and if we assume = 1 it is 
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m 

= _E 
C A P 



(9) 


What Eq. (9) tells us is that droplets splashed back from the surface of an 
airfoil will have the momentum defect that they bring into the air flow about 
the airfoil adjusted within a distance of the order of a from the surface. 

We may look at this length then as a momentum-defect deposition or adjustment 
length. We may express this length in terms of the splash droplet diameter, 
d p , by noting that for a sphere 



§ p w d p 


(10) 


where P w is the density of water. Putting Eq. (10) into Eq. (9) we obtain 


Using the diameter d p given in (6) we find finally that the momentum-defect 
adjustment length is 


l 
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16 °w 

3f p V 2 
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( 12 ) 


What will be important in considering the effect of the momentum defect 
splashed back by the impinging droplets will be the ratio of this length to a 
typical dimension of the wing or airfoil in question. We therefore expect 
that an appropriate parameter for scaling the effect of heavy rain will be 
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where c is the chord of the airfoil. We can get rid of the constant 3f/16 
and take as a parameter simply 
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P a V’ c 
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(14) 


This parameter (the deposition length parameter) is a mixed Weber number con- 
taining the density of air and the surface tension of water. 


It is useful in order to understand the relative importance of rainfall 

rate and N on C , to consider the sketches in Figure 2.2a and 2.2b. In 
D L 

max 

Figure 2.2a we show for a fixed N D (thus a fixed momentum-adjustment length) 
an extremely simplified picture of what increasing rainfall rates will do to 
the velocity distribution on the upper surface of an airfoil at angle of at- 
tack. Clearly as the rainfall rate is increased, a larger momentum-defect 
will be realized, and with increasing defect will come less and less ability 
of the boundary layer on the upper surface to recover pressure. Thus, we 
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.2a Effect of increasing rainfall rate on 

near-surface velocity distribution at a 
fixed value of Np, i.e. fixed momentum- 

defect adjustment length relative to chord. 
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might expect the of an airfoil to fall off (perhaps linearly) as rain- 

max 

fall rate is increased at a fixed value of N^. 

On the other hand, if we consider what happens for a fixed rainfall rate 
as is increased, we must see something like what is depicted in Figure 
2.2b. In this case, as is increased, the momentum adjustment length be- 
comes smaller and smaller; thus a given momentum deficit that is associated 
with the rainfall rate must be carried in a thinner and thinner layer. The 
result, at least as far as the flow at the surface is concerned, is equivalent 
to increasing the rainfall rate. We would expect to find then that the effect 
of increasing the chord of an airfoil at a fixed rainfall rate might be 
similar to the effect of increasing rainfall rate for a fixed chord. 
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Figure 2.2b Effect of increasing Np on near surface 

velocity distribution at fixed rainfall rate. 
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3. Initial Examination of Data 


-# 
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A small amount of data on C as a function of rainfall rate and N_, was 

Li U 

max 

given to the authors by Mr. Earl Dunham of NASA/Langley. These data are shown 

in Figure 3.1. A fair amount of liberty has been taken in fairing curves 

through the data on the basis of our feeling that the effect of rainfall rate 

on C L might be linear. The fairing was also biased by our notion that the 
max 

effects of Np and rainfall rate might be similar. 

With great trepidation the results of the initial fairing of these data 
are shown in Figure 3.2, as a function of with rainfall rate as a 
parameter. The C T data are plotted in terms of 1-C T / ( C T ) rt , where 

Li Li Li U 

max max max 

(C^ )q is the value of for zero rainfall rate, 

max max 

On this figure we have also indicated the range of N D that will be found 
for large jet transport aircraft during landing operations and the range of N Q 
for which water impingement tests have been carried out at MIT. 

First of all, it should be noted that Reynolds number is a parameter that 

must be considered, perhaps by a consideration of the ratio of i to the 

c 

boundary layer thickness. However, it is now believed best to use as given 
here as the essential parameter and consider Reynolds number effects 
separately. 

In view of the possibility that effects are similar to rainfall rate 
effects as pointed out above, an initial conjectural extrapolation of the data 
to the Np values that might be associated with large scale transports was made 
and is given in Figure 3.2 for a rainfall rate of 75 mm/hr (3 in/hr). This 
extrapolation indicated that such a rate might cause a 15 to 20 percent reduc- 
tion of C L . Clearly the implication of this hypothetical scaling poses a 
max 
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Figure 3.1 NASA data on C L as a function of rainfall rate with N as 

max D 

parameter. 
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Figure 3.2 Conjectural plotting of NASA data on 
with rainfall rate as a parameter. 






sufficiently severe operational problem so that a large scale wing should be 
tested in a rain environment. 

In addition it was clear than an attempt to calculate the effect of 
droplet splash-back on the maximum lift of an airfoil would be a most useful 
exercise. In what follows we will describe the program that was used to make 
such calculations and the results of these calculations. 


4. Program Development and Use 


The calculation of the maximum lift of a two-dimensional airfoil even 
without rain is not yet a science. Therefore, rather than try to calculate 
the effect of rain on maximum lift, the effect on the position of separation 
for an airfoil near maximum lift was studied. An NACA 64-210 airfoil at an 
angle of attack, a, of 12° was chosen. 

To determine the point of separation, a program called ARB, which has 
been in use at A.R.A.P. for several years (6), was used. ARB, which employs 
full second-order closure models for turbulence, computes the boundary layer 
(using standard boundary-layer assumptions) on an axisymmetric body of ar- 
bitrary shape, whether rotating or not, as well as on two-dimensional bodies. 
It handles compressible flows, solving for the extra second-order correlations 
that exist because of fluctuations in temperature and density, as well as 
constant-density flows. Since the present study concerns low-speed flows, the 
constant density mode was used. 

A special version of ARB was developed for this project in order to cal- 
culate the effect of particles (water droplets) as they splatter back into the 
boundary layer after hitting the surface of an airfoil. A description of the 
analysis underlying the modification of ARB follows. 

4.1 Motion of a particle in a boundary layer 

In analyzing the motion of a particle in a boundary layer, it is assumed 
that the variation of the mean flow with x, the streamwise direction, is neg- 
ligible. It is also assumed that the flow is two dimensional and that the 
effect of particle motion in the third dimension can be neglected after 
averaging over a large number of particles. 


Consider then a particle moving with velocity v = (u ,v ) in a mean flow 
-*■ P P P 

v = (u,v) (functions of y only). In terms of the relative velocity, 

— ► — ► -+• 

v = v-v , the force on the particle is 
r p 

~ C A p Iv |v 
2 D p a 1 r 1 r 

where C D is the drag coefficient, A p is the cross-sectional area of the par- 
ticle, and p is the density of air. Then the acceleration of the particle is 
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VA -D 

Tt ' ""2m 'W ' T~ l v r' v r 
P c 


where m^ is the mass of the particle and 


c A p 
P a 


To evaluate the drag coefficient, the particle is assumed to be a sphere 

of diameter d . In terms of the Reynolds number 
P 


| v |d 
r t 

R ■= — 1 1 

p V 


where v is the kinematic viscosity of air, may be approximated by 
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(See Figure 1.5 in (7).) For a sphere, 
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where p 


w 


is the density of water, and so 


C D = 


32v 


p w 

— + 1 
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The acceleration can be rewritten 



dt 
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K p 



= D„(v-v ) 
K p 


where 
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In order to evaluate the effect of the particles on the air flow it is neces- 
sary to know Vp as a function of y. Since 


v 

P 


= & 
dt 


the acceleration equations can be written 
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dy v p p' 


dv D 

— = — ( v -v ) 
dy v p 
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Thus if the initial values of u and v as a particle comes off the surface 

P P 

are known, and if the flow conditions are known as a function of y, these 
equations can be integrated and the force evaluated. 

If there are n^ particles per unit volume so that the liquid water con- 
tent is p * n w , the x component of the total force exerted by the particles 
P P P 

on the air, equal and opposite to the drag of the particles, is 

p pW u) 

per unit volume. This term was added on the right side of the mean momentum 
equation (x component) in ARB. No way was found to incorporate the effect of 
the vertical component of the force from the particles on the fluid in the 
context of the boundary-layer assumptions on which ARB is based, since these 
assumptions imply that the vertical component of the mean momentum equation is 
automatically in balance and the vertical component of velocity is computed 
from the continuity equation. 

However it was recognized that if £ is very small, a particle will give 

0 

up its momentum to the air almost immediately so that the effect can be simu- 
lated by applying a blowing boundary condition with the velocity at the wall, 
v , given by 


v 

a w 


= p v (0) 
P P 


The results of these runs are shown in Section 5 as limits for £ /c equal 

cm 

zero although they were actually run with £ /c of the order of 2xl0 -5 . 

cm 


The rate at which the wake of a particle receives energy is 


d 'V 2 

dt (m p“F" ) = 


dv C_m 

- r Dp,-*, 
- m v • ~rr~ = : v 

p r dt d 1 r 
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Therefore the rate at which the fluid receives energy per unit volume is 


p p d k 'V’ 

It is assumed that a fraction, e, of this energy produces turbulence, the 
rest going into heat. Therefore a source term equal to 


f p p °K 'V 2 


was included in the equations for u'u' , v'v' , and w'w' which together con- 
stitute the turbulent energy. For all the runs described in Section 5, e was 
taken to be 1/2. 

4.2 The variation with streamwise distance 

The analysis above assumes that the variation with x is negligible; 
however as the boundary-layer calculation moves downstream from the stagnation 
point the flow velocity v(y) changes (partly as a result of the force calcu- 
lated above) and the characteristics of the particles splashed back also 
change . 

Let v q be the normal component of the velocity of a raindrop as it im- 
pinges on the airfoil surface. Then, neglecting the settling speed of the 
raindrops in comparison with the aircraft speed, u , 


v = u sin 0 
o «® 


where 0 is the local angle of the airfoil surface. 


19 




The particles are assumed to leave the surface normally at a given proportion 
of v , say p. Then 

u p (0) = 0 

v (0) = pv = pu sin0 = v sine 
p o r «> pm 

where v pm = pu^. It is assumed that no particles are splashed back where 0 is 
zero or negative. 

It is shown in Section 2 that 2 is inversely proportional to the square 
of v q so it can be written 


2 
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2 /sin 2 0 
cm 
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where 2 is the minimum value of 2 . 
cm c 
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4.3 Pressure distribution 


ARB requires the free-stream pressure distribution or the free-stream 
velocity distribution (from which the pressure distribution can be calculated) 
as input. For the runs reported in the next section the inviscid flow over 
the airfoil was calculated by the method developed by J. L. Hess. The results 
of these calculations were found to be a little noisy so they were passed 
through a smoother before being supplied to ARB as outer boundary conditions. 

The distribution of u , the free-stream velocity, as a function of x, 

e 

measured along the surface from the stagnation point, is shown in Figure 4.1. 
Figure 4.2 is a diagram of the airfoil at a - 12° with the stagnation point, 
the leading edge (defined here as the point where 0 = 90°), the point of mini- 
mum pressure (maximum u g ), and the point of zero slope (0 = 0) indicated. It 
will be noted that the point of minimum pressure is almost at the leading edge 
and that on the upper surface a region of large adverse pressure gradient oc- 
curs between the leading edge and the point of zero slope. In this region we 
might expect a very heavy rain to cause the flow to separate and hence cause 
premature stall of the airfoil. 





5. Computations and Results 


Many calculations were made of the behavior of the boundary layer for a 

number of conditions of rainfall rate (represented by p ) , particle size 

P 

(represented by H ), initial particle velocity (represented by v ), and 
cm pm 

Reynolds number. Figure 5.1 shows typical behaviors of calculated skin fric- 
tion coefficient, c., for three values of p /p Q with H /c = 0.001, v /u = 

f p a cm pm 00 

0.4, and a Reynolds number of 3.5xl0 7 . This is a Reynolds number typical of a 

747 aircraft near landing speed. The curve for p /p =0 shows that, as ex- 

P a 

pected, the aircraft does not stall and the point of separation (c^ = 0) is 

near the trailing edge of the airfoil. The curve for p /p =0.1, repre- 

P a 

senting a high rainfall rate, shows that the separation point has moved 

forward somewhat so we might expect some drop-off in C L but the airfoil has by 

no means stalled. Note however that the dip in c f in the region of strong 

adverse pressure gradient has become more pronounced. The third curve, for a 

still higher rate, p /p =0.12, shows that the separation point has jumped 

p a 

forward to that region and we would expect that the airfoil has stalled. 

Figure 5.2 shows the value of x/c at which separation occurs (plotted 

horizontally) as a function of p /p (plotted vertically) for the three runs 

p a 

shown in Figure 5.1 and other runs all at the same i /c, v /u , and Reynolds 

cm pm 00 

number. It is clear from this figure that the critical value of p /p , that 

pa 

required to cause early separation (which we assume is equivalent to stall), 
is 0.113. 

Repeating the process described above for different values of a cm / c * 

still at the same initial velocity and Reynolds number, we obtain the results 

shown in Figure 5.3. The open circles are the calculated values of p /p re- 

p a 

quired to trigger early separation at various values of H^/c. When results 
such as these were first encountered they caused some consternation, for cer- 
tainly we did not expect the rainfall rate required for early separation to 
increase as the momentum defect adjustment length, H c , decreases. The ques- 
tionable points are connected by a dashed curve in the figure. An analysis of 


004 



25 


Figure 5.1 Variation of the local coefficient of skin friction as a 

function of distance from the stagnation point, with rainfall 
rate (represented by p^) as a parameter, on an NACA 64-210 

airfoil at a = 12° for l /c = 0.001, v /u = 0.4, and 
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Figure 5.3 Variation of the critical value of p p as a function of t^/c for 

an NACA 64-210 airfoil at a = 12° with v /u =04 and 
„ pm 00 

Re = 3.5xl0 7 . 




this phenomenon showed it to be related to the fact that as 2 becomes small 

0 

compared to the boundary- layer thickness, the effect of the particles on the 
vertical momentum in the boundary layer needs to be included in the calcula- 
tions. As mentioned in Section 4.1, this is not possible in the ARB code, 
but a way was found to get around this problem, and it certainly is a problem, 

by using the surface gas injection feature of ARB for the limit 2. =0. The 

c 

result of such a calculation is shown by the filled circle for 2 /c = 0. For 

cm 

the purposes of estimation, the authors, for want of any better method, have 
faired a smooth curve shown by the dot-dash curve in Figure 5.3, between this 
point and the other points we believe to be valid. 

Curves such as that shown in Figure 5.3 have been calculated for a number 

of Reynolds numbers and two values of the ratio v /u . Summary plots of this 

pm 00 

nature obtained from many runs are shown in Figures 5.4 and 5.5. 


Additional scales have been included on these plots. 
14 of Section 2 are combined, we find 


If equations 13 and 


N = — — — 

D 3f 2 

cm 


where 2 has been used instead of 2 because the value of V meant in the 
cm c 

definition of N„ (= p V 2 c /O ) is here taken to be the speed of the aircraft, 
Daw 

u . (See Section 4.2.) Taking the value of f, the portion of the kinetic 
00 

energy of the drops impinging on the surface converted to surface energy, to 
be 1/2, we can write 


N S-™_ 
D 2 /c 
cm 


Values of thus calculated are indicated as additional horizontal scales in 
Figures 5.4 and 5.5. The additional vertical scales are rainfall rates, ob- 
tained from p /p under the following assumptions. The settling velocity of 
P ® 

the raindrops is taken to be about 5 m/sec and it is assumed that half the 
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Figure 5.5 Variation of the critical value of p p or rainfall rate as a 

function of t /c or N , with Reynolds number as a parameter 
cm u 

for an NACA 64-210 airfoil at a = 12° with v /u w =0.6. 


water impinging on the wing is splashed back. Then it can be shown that at 

standard conditions the rainfall rate in mm/hr is approximately 4xl0 4 p /p . 

pa 

Examining Figure 5.4, for v /u * 0.4, we note that, as expected, for a 

pm 00 

low Reynolds number (6.5xl0 5 ) the airfoil stalls even when there is no rain- 
fall. As the Reynolds number is increased to the level of the NASA tests 
mentioned above (2.4xl0 6 ), the rainfall required to cause early separation 

decreases as l /c decreases, i.e. as N_ increases. For Reynolds numbers 
cm D 

typical of transport aircraft (1.8xl0 7 and 3.5xl0 7 ) the calculations indicate 
that the increased Reynolds number offsets the effect of pure size on so 
that a somewhat larger rainfall rate is required to induce early stall for a 
747 aircraft than was found for conditions of the NASA tests. This is an in- 
teresting result and is at variance with the pessimistic results of the 
authors' first extrapolation of the NASA data (see Figure 3.2) 

Figure 5.5 shows similar results for = 0.6. If was anticipated 

that this parameter, which relates the normal velocity of the incoming drop to 
the initial velocity of the droplets resulting from its fracture, would be a 
very sensitive parameter in the interaction between the droplets and the 
boundary-layer. A comparison of Figures 5.4 and 5.5 shows this to be the 
case. For example, the rainfall rate required to induce early separation for 
the conditions of the NASA tests dropped from about 1400 mm/hr to about 800 


6. Discussion and Recommendations 


The computations that have just been presented are believed to give the 
general trends associated with the phenomenon of premature stall due to cou- 
pling of the splash-back droplets with the boundary layer. Clearly the 
numerical values of rainfall rate to produce premature separation of a given 
Reynolds number and mixed Weber number (the deposition number Np ) are not 
exact. However, we believe this to be of the right order of magnitude. 

If this conclusion is correct it follows that for rainfall rates ap- 
proaching 500 mm/hour or greater, commercial aircraft might be subject to rain 
induced premature stall. The calculations also indicate that both Reynolds 
number and scale are important in pinning down the values of rainfall rate 
that will result in early separation due to ejecta coupling with the boundary 
layer. In view of these facts and in view of the fact that other extremely 
complex phenomena also have an effect on boundary layer behaviors in the 
presence of heavy rain, it is clearly necessary to conduct full scale tests if 
a quantitative assessment of heavy rain hazard is to be made. This is true 
for simple wings and is certainly even more true of the flapped and slotted 
wings typical of commercial aircraft. 

In regard to analysis, although it is believed that certain trends that 
have been calculated with the extended ARB program are correct, it has been 
shown that, if more exact computations are required, a method more powerful 
than conventional boundary layer theory must be applied. Since such analysis 
will be expensive, it seems desireable to hold off any further development of 
analytical capability until such time as both full scale and model tests of 
the same wing configuration can be completed so that the importance of heavy 
rain effects on wing performance can be evaluated more precisely. 


References 


1. Rhode, R. V., "Some Effects of Rainfall on Flight of Airplanes and on 
Instrument Indications," NACA Tech. Note 803, 1941. 

2. Haines, R. A. and Luers, J. K. , "Aerodynamic Penalties of Heavy Rain on 
a Landing Aircraft," J. Aircraft, Vol. 20, No. 2, February 1983. 

3. Dunham, E. R., Bezos, G. M., Gentry, C. L. and Melson, E., "Two- 
Dimensional Wind Tunnel Tests of a Transport-Type Airfoil in a Water 
Spray," AIAA-85-0258 . 

4. Donaldson, C. duP., "On the Scaling of the Effect of Heavy Rain on the 
Maximum Lift of Airfoils, A.R.A.P. Technical Memorandum 85-33, August, 
1985. 

5. Bilanin, A. J., "Scaling Lavs for Testing of High Lift Airfoils Under 
Heavy Rainfall," presented at AIAA 23rd Aerospace Sciences Mtg. , January 
1985, Reno, Nevada, AIAA-85-0257. 

6. Sullivan, R. D. , and Varma, A. K. , "ARB: A Program to Compute the 

Turbulent Boundary Layer on an Arbitrary Body of Resolution," A.R.A.P. 
Report No. 317, January 1978. Also: Sullivan, R. D. , "ARB: A 

Supplementary Manual," A.R.A.P. Tech Memo 80-16, September 1980. 

7. Schlichting, H., Boundary Layer Theory, Fourth Edition, McGraw-Hill Book 
Co., Inc., translated by J. Kestin, 1960. 



Standard Bibliographic Page 


1. Report No. 

NASA CR-178248 


4. Title and Subtitle 


2. Government Accession No. 


3. Recipient’s Catalog No. 


5. Report Date 

March 1987 


6. Performing Organization Code 


8. Performing Organization Report No. 

ARAP Report No. 597 


10. Work Unit No. 


11. Contract or Grant No. 

NAS1-18088 


13. Type of Report and Period Covered 

Contractor Report 


14. Sponsoring Agency Code 

505-45-13-01 


The Effect of Heavy Rain on an Airfoil at High Lift 


7. Author(s) 

Coleman duP. Donaldson and Roger D. Sullivan 


9. Performing Organization Name and Address 

A.R.A.P. Group 
Titan Systems, Inc. 
Princeton, New Jersey 08540 


12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 


15. Supplementary Notes 

Langley Technical Monitor: R. Earl Dunham, Jr. 


16. Abstract 

Although the effect of heavy rain on aircraft performance was discussed 
as early as 1941 by Rhode (1), no serious studies of the relationship of heavy 
rain to aircraft safety were made until Luers suggested in 1981 that the 
torrential rain that often occurs at the time of severe wind shear might 
substantially increase the danger to aircraft operating at slow speeds and 
high lift in the vicinity of airports. While Luer's ideas were not published 
until early 1983 (2), appropriate measures were taken by NASA to study the 
effect of heavy rain on the lift of wings typical of commercial aircraft. 

These tests, reported by Dunham, Bezos, Gentry, and Melson (3)» were the 
subject of a number of discussions between the senior author of this report 
and Mr. Earl Dunham of NASA during the fall of 1984. One of the aspects of 
these tests that seemed confirmed by the data was the existence of a "velocity 
effect" on the lift data. The data seemed to indicate that when all the 
normal non-dimensional aerodynamic parameters were used to sort out the data, 
the effect of velocity was not accounted for, as it usually is, by the effect 
of dynamic pressure. Indeed, the measured lift coefficients at high lift 
indicated a drop-off in lift coefficient for the same free-stream water 
content as velocity was increased. 
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